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ABSTRACT 

HD 89744 is an F7 V star with mass 1.4 M , effective temperature 6166 K, age 2.0 Gy 
and metallicity [Fe/H]= 0.18. The radial velocity of the star has been monitored with the AFOE 
spectrograph at the Whipple Observatory since 1996, and evidence has been found for a low mass 
companion. The data were complemented by additional data from the Hamilton spectrograph 
at Lick Observatory during the companion's periastron passage in fall 1999. As a result, we have 
determined the star's orbital wobble to have period P — 256d, orbital amplitude K = 257 m/s, 
and eccentricity e = 0.7. From the stellar mass we infer that the companion has minimum mass 
77i2 sin i = 7.2 Mjup in an orbit with semi-major axis — 0.88 AU. The eccentricity of the orbit, 
among the highest known for extra-solar planets, continues the trend that extra-solar planets 
with semi-major axes greater than about 0.15 AU tend to have much higher eccentricities than 
are found in our solar system. The high metallicity of the parent star reinforces the trend that 
parent stars of extra-solar planets tend to have high metallicity. 



Subject headings: planetary systems 
techniques:radial velocities 

1. INTRODUCTION 



stars:low-mass, brown dwarfs — stars: individual (HD 89744) 



We report on the detection of a massive 
(rr^sini = 7.2 Afjup) planet in a highly ellip- 
tical (e = 0.7), 256 day orbit about the star 
HD 89744 (HR4067, HIP 50786), from radial ve- 
locity variations which reveal Keplerian motions 
of the star. Observations were carried out from 
1996 through 1999 using the Advanced Fiber Op- 
tic Echelle (AFOE) spectrograph (Brown et al. 
1994; Nisenson et al. 1998), a bench-top spectro- 
graph located at the Whipple Observatory 1.5m 
telescope, and also with the Hamilton spectro- 
graph at the Lick Observatory CAT and Shane 
telescopes, in November and December of 1999. 

The AFOE spectrograph is designed primarily 
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for precise radial velocity studies of the seismol- 
ogy of bright stars, and of reflex motions of stars 
due to planetary companions. Long term stability 
of the velocity reference is provided by use of an 
iodine (J 2 ) cell (Butler et al. 1996). The AFOE 
determines radial velocity variations induced by 
planetary companions with a precision and long- 
term accuracy of approximately 10 m/s. On the 
order of 100 relatively bright stars (m v < 7) have 
been monitored for this purpose since 1995. 

Since 1995 when the planetary candidate oribit- 
ing the star 51 Pegasus was detected (Mayor & 
Queloz 1995), some 29 additional candidates have 
been detected by several groups, all from Doppler 
shifts measured using precise radial velocity tech- 
niques (Marcy & Butler 1998; Mayor et al. 1998; 
Noyes et al. 1997; Cochran et al. 1997). 

HD 89744 (F7 V) was added to the AFOE ob- 
serving list in early 1996, based on its relatively 
low chromospheric emission as measured with the 
Mt. Wilson "HK" chromospheric activity moni- 
toring program (Baliunas et al. 1995). AFOE 



observations have been obtained regularly since 
then, and indicated the presence of a planet with 
a highly eccentric orbit. However data near the 
companion's periastron, critical to an accurate de- 
termination of the orbital parameters, were not 
obtained until late 1999. Between October and 
December 1999, while the companion was near pe- 
riastron, observations were made at Lick Observa- 
tory as well as with the AFOE, to ensure good 
phase coverage. The data points taken with the 
Lick CAT and Shane telescopes agree extremely 
well with the AFOE data, and thus provide a con- 
firmation of the detection along with a precise de- 
termination of the ellipticity of the planet's orbit. 

2. PROPERTIES OF THE HOST STAR, 
HD 89744 

HD 89744 is an F7V star at a Hipparcos- 
determined distance of 39.0 parsec. It is listed 
as a constant star in the Hipparcos catalog . The 
star has absolute magnitude M v — 2.78 and color 
B — V = 0.531 (Perryman 1997). Comparing 
its position in the color-magnitude diagram with 
predictions of stellar evolution calculations, Pri- 
eto & Lambert (2000) determine it to have mass 
M = (1.34±O.O9)M , radius R = (2.14±O.1)-R , 
and effective temperature T c ff= (6166 ± 145) If . 
Independently, Ng & Bertelli (1998) determine 
its mass to be M = (1.47 ± 0.01)M©. For the 
purposes of this paper we adopt the average of 
the two masses listed and an uncertainty given by 
their spread: M = (1.4 ± 0.09)M Q . 

The metallicity of HD 89744 has been deter- 
mined to be [Fe/H]= 0.18 by Edvardsson et al. 
(1993). Its age is determined by Ng & Bertelli 
(1998) to be 2.04 ±0.10 Gy. 

From rotational modulation of the Ca II flux, 
Baliunas, Sokoloff, & Soon (1996) determined that 
the rotation period of the star is P ro t = 9 days. 
This, combined with the above-mentioned radius 
of the star, implies an equatorial velocity v = 12 
km/sec. Rotational broadening of the spectrum 
implies vsini = 8 km/s (Bcrnacca & Perinotto 
1970; Uesugi & Fukuda 1970). This further im- 
plies that the star's rotational equator is inclined 
by about 40° to the plane of the sky. 

Speckle observations (McAlister et al. 1989) 
have not revealed any indication of a stellar com- 
panion to HD 89744. 



Table 1 summarizes the relevant parameters of 
the star HD 89744. 

Table 1: Parameters for HD 89744. 



Parameter Value 



Mass [M e ] 


1.4 


± 


0.09 


Radius [Rq] 


2.14 


± 


0.1 


T off [K] 


6166 


± 


145 


M v 


2.78 






Age [Gy] 


2.04 






[Fe/H] 


0.18 






P mt [days] 


9 






vsini [km/s] 


8 







3. OBSERVATIONS 

3.1. Instrumentation and Data Reduction 

The AFOE is a bench-mounted, fiber-fed, high- 
resolution cross-dispersed echelle spectrograph de- 
signed for high radial velocity precision and sta- 
bility both on short time scales (better than 1 
m/s over hours, for asteroseismology) and long 
term (approximately 10 m/s, for radial velocity 
exo-planet searches, through use of an iodine ab- 
sorption cell). The AFOE is located at the 1.5m 
telescope of the Whipple Observatory on Mt. Hop- 
kins, Arizona. A more complete description of 
the AFOE is available in Brown et al. (1994) and 
Nisenson et al. (1998). The AFOE exo-planet sur- 
vey program monitors the radial velocity of about 
100 stars brighter than m v — 7.5, with an accu- 
racy of 10 - 15 m/s for integrations with a signal- 
to-noise ratio of 100 to 150. Most observations 
consist of three consecutive exposures, primarily 
to limit cosmic ray contamination by keeping the 
exposure times short. 

Our data reduction methodology is conceptu- 
ally similar to that described by Butler et al. 
(1996), but differs in details. Echelle images are 
dark subtracted and one dimensional spectra ex- 
tracted and corrected for scattered light, then flat- 
fielded using the spectrum of a tungsten lamp. For 
each of the six spectral orders that contain strong 
iodine lines a model is adjusted to match the 
observed star-plus-iodine spectrum in the least- 
squares sense. This model is computed using a 
Dopplcr-shiftcd high SNR spectrum of the star 
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alone, plus a very high resolution high SNR spec- 
trum of the iodine cell. The model incorporates 
the sought-for Doppler-shift of the star as well as 
mechanical drifts within the spectrograph, the in- 
strumental wavelength solution, the instrumental 
resolution profile, and a residual scattered light 
correction. The resulting radial velocities, after 
correction for the motion of the telescope relative 
to the solar system barycenter, are averaged for all 
six orders and the three consecutive observations. 
The scatter around the mean provides an estimate 
of the uncertainty. The root-mean-square (RMS) 
velocity of several radial velocity standard stars is 
commensurate with this estimate of uncertainty. 
The Doppler analysis for the Lick observations is 
described in detail in Butler et al. (1996). 

3.2. Observations and Orbital Fit 

AFOE observations of HD 89744 were obtained 
on 74 separate nights between December 1996 and 
December 1999. Additional observations on 14 
nights during November and December 1999 were 
also taken at the Lick Observatory CAT and Shane 
telescopes, to ensure complete phase coverage near 
the November 1999 periastron passage. The data 
with their uncertainties are plotted in Figure 1. 

The RMS of the zero-averaged velocities for the 
AFOE data alone is 129.1 m/s, while the averaged 
uncertainty on these observations is 15.2 m/s; the 
corresponding reduced x is 10; where the reduced 
X is given by 

where Oi are the observations, at the observational 
uncertainties, Mj the corresponding values of the 
model, N the number of observations and N, m 
the number of adjustable parameters in the model. 
After fitting a Keplcrian orbit in the least-squares 
sense, the RMS of the residuals is 26.2 m/s, corre- 
sponding to a reduced x OI 1-7. If we include the 
Lick observations, by adding one additional free 
parameter to account for the arbitrary offset be- 
tween the two data sets, the RMS of the residuals 
is 20.5 m/s, or a reduced x of 1.6. (One velocity 
point was rejected based on a 3-pass 3-cr rejection 
algorithm.) The resulting orbital parameters are 
given in Table 2. These are almost identical to the 
parameters obtained when using only the AFOE 



observations, except that the combined data set 
leads to smaller uncertainties, primarily on the 
amplitude but also on the eccentricity. The orbital 
fit phase plot and the residuals to the fit are shown 
in Figures 2 and 3. The periodic variation of the 
radial velocity, and the highly eccentric character 
of the orbit, are evident. For a stellar mass of 1.4 
Mq, the orbital elements imply a minimum mass 
for the companion of m,2 sini = 7.2 Mjup, and a 
semi-major axis of 0.88 AU. 

Table 2: Orbital parameters resulting from a Ke- 
plcrian fit to the combined data set. The last two 
rows of the table use the adopted value mi = 1.4 

M Q 



Parameter 


Value 




P [day] 


256.0 


± 


0.7 


K [m/s] 


257 


± 


14 


e 


0.70 


± 


0.02 




195 


± 


3 


T [JD-2,450,000] 


994 


± 


2 


f(mi,m 2 ,i) [M Q ] 


1.64 10" Y 






ai ami [AU] 


0.0043 






RMS (residuals) [m/s] 
reduced x 


20.5 
1.6 






iV To t, iV Re j 


88, 




1 


a 2 [AU] 

77i2 sin i [Mjup] 


0.88 
7.2 







4. DISCUSSION 

The radial velocity data shown in Figure 2 are 
unambiguous in revealing a periodic radial veloc- 
ity variation of HD 89744, which can be fit well 
by a Keplerian orbit. Since observations obtained 
with two different radial velocity instruments at 
two different telescopes yield exactly the same ra- 
dial velocity variations within their respective in- 
strumental errors, the evidence is compelling that 
the measured velocities are real variations on the 
star. There is no known way a stellar signal in a 
late F-type main sequence star could mimic a Ke- 
plerian orbital signature with such a long period 
and large amplitude, and with such a large eccen- 
tricity. Hence we are driven to the interpretation 
that the star is orbited by a low-mass companion, 
HD 89744 b, (m 2 sin7 = 7.2M Q ), in an orbit with 
semi-major axis 0.88 AU and eccentricity 0.70. 
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Fig. 1. — Radial velocity observations of HD 89744 with their respective uncertainties. Filled circles denote 
AFOE observations, and open circles Lick observations after applying an offset determined by the combined 
orbital fit. The Keplerian orbital fit is indicated by the dotted line. 



The residuals to the orbital fit are larger than 
would be expected from internal errors in the data; 
this is true both for the AFOE data and the Hamil- 
ton echelle data. However, Saar et al. (1998) con- 
clude that for a typical F stars with this rotation 
period, the velocity jitter induced by stellar mag- 
netic activity and inhomogeneous convection is ap- 
proximately 10 m/s. Adding this jitter in quadra- 
ture to the uncertainties lower the reduced \ to 
1.2. The residuals display a long term trend (« 15 
m/s/year) that is marginally significant. While it 
might be caused by a residual instrumental drift in 
the AFOE data, we can not rule out that it might 
be due to a distant companion; further observa- 
tions over a longer baseline are required. 

The orbital eccentricity of the companion to 
HD 89744 is among the highest planetary ec- 
centricities known. Only two other planets, 
16CygB b (e = 0.68, a 2 = 1.70 AU; Cochran 
et al. 1997) and HD 222582 b (e = 0.71, a 2 = 1.35 
AU; Vogt et al. 1999) have comparable eccentric- 



ities. HD 89744 b, with a 2 = 0.88 AU, has the 
smallest semi-major axis of the three. At perias- 
tron it dips to within 0.26 AU, still well outside 
a periastron distance which could lead to tidal 
circularization within the stellar lifetime. 

The discovery of the first highly eccentric 
planet, 16CygB b, led to the suggestion (Mazeh 
et al. 1997; Holman et al. 1997) that its eccentric- 
ity may have been "pumped up" by the influence 
of a nearby companion star, 16CygA. However, 
neither HD 89744 nor HD 222582 is orbited by a 
stellar companion. Thus a different explanation is 
required for their high eccentricities, an explana- 
tion that might also apply to 16CygB b. 

The high orbital eccentricity of the HD 89744 
system continues the trend that planetary-mass 
companions whose semi-major axes are greater 
than about 0.15 AU tend to have a broad range 
of eccentricities, with no apparent trends of ec- 
centricity with mass or semi-major axis. This cir- 
cumstance must be explained by any successful 
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Fig. 2.— Phase plot of the Keplerian orbital fit. Filled circles denote AFOE observations, open circles Lick 
observations. 



planetary formation and migration scenario. As 
noted by others (e.g., Vogt et al. 1999; Marcy et 
al. 1999), this causes difficulties with a number of 
proposed mechanisms for planetary formation and 
migration. 

The values of vsmi, P rot , and R for HD 89744 
given in Table 1 imply an inclination of the stellar 
rotational equator of i = 42° (i.e., sini = 0.66). 
Moreover, if we assume that the orbit is copla- 
nar with the star's equatorial plane, we infer that 
TO2 = 10.8 Mjup- The astrometric orbital ampli- 
tude would be 0.17 mas, too small for Hipparcos 
detection but within the range of next-generation 
astrometric missions. 

The mass of 10 Mjup suggested by the stel- 
lar rotational data is near the upper limit of 
masses associated with extra-solar giant planets 
(e.g., Marcy et al. 2000). If such large masses hold 
up to further investigations, then theoretical un- 
derstanding of the origin and evolution of extra- 
solar giant planets must be able to accommodate 
a mass range spanning at least values between 0.5 



Mjup and 10 M JUP . 

The metallicity of HD 89744, [Fe/H]= 0.18, is 
substantially higher than the mean for nearby 
sun-like stars (Favata et al. 1997; Gonzalez 1998). 
HD 89744 was placed on the AFOE observing list 
without reference to its metallicity; therefore the 
association of its high metallicity with the pres- 
ence of a planet is not a selection effect. This as- 
sociation continues the trend, already noted (e.g., 
Gonzalez & Law 1999, and references therein) 
and references therein that the metallicity of stars 
with planets tends to be higher than that of stars 
without planets. 

We are grateful to the Mt. Hopkins observ- 
ing and support staff, especially Ed Bennet, 
Perry Berlind, Mike Calkins, Ted Groener, Robert 
Hutchins, Jim Peters and Wayne Peters; we also 
acknowledge observing with the AFOE by Scott 
Horner and Ted Kcnnelly. We are very grateful 
to Adam Contos for his dedicated efforts in ob- 
serving and reducing a significant fraction of data 
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Fig. 3. — Residuals to the Kcplerian orbital fit. Filled circles denote AFOE observations, open circles Lick 
observations. The dotted lines illustrate the 3-sigma rejection bounds. 
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